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Abstract Because of their deep gravitational potential wells, clusters of galaxies retain all 
the metals produced by the stellar populations of the member galaxies. Most of these met- 
als reside in the hot plasma which dominates the baryon content of clusters. This makes 
them excellent laboratories for the study of the nucleosynthesis and chemical enrichment 
history of the Universe. Here we review the history, current possibilities and limitations of 
the abundance studies, and the present observational status of X-ray measurements of the 
chemical composition of the intra-cluster medium. We summarise the latest progress in us- 
ing the abundance patterns in clusters to put constraints on theoretical models of supernovae 
and we show how cluster abundances provide new insights into the star-formation history of 
the Universe. 
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1 Introduction 

Clusters of galaxies are excellent astrophysical laboratories, which allow us to study the 
chemical enrichment history of the Universe. They have the deepest known gravitational 
potential wells which keep the metals produced in the stellar populations of the member 
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galaxies within the clusters. About 70-90 % of the baryoni c mass content of clu sters of 
galaxies is in the form of hot (10^-10* K) X-ray emitting gas tettori & FabiarJ 19991) . In this 
hot intra-cluster medium (ICM) the dominant fraction of cluster metals resides. To the extent 
that the stellar populations where the cluster metals were synthesised can be considered 
representative, the metal abundances in the ICM provide constraints on nucleosynthesis and 
on the star formation history of the Universe. 

It is remarkable that all the abundant elements, that were synthesised in stars after the 
primordial nucleosynthesis, have the energies of their K- and L-shell transitions in the spec- 
tral band accessible to modem X-ray telescopes. Most of the observed emission lines in 
the ICM arise from the well understood hydrogen and helium like ions and their equivalent 
widths can be, under the reasonable assumption of collisional equilibrium, directly con- 
verted into the elemental abundance of the corresponding element. Complications arising 
from optical depth effects, depletion into dust grains, extinction, non-equilibrium ionisation 
are minimal or absent. Therefore, abundance determinations of the hot ICM are more robust 
than those of stellar systems, H ii regions or planetary nebulae. This relatively uncomplicated 
physical environment makes the ICM an attractive tool for studies o f the chemical enrich- 
ment (for the theoretical progress in metal enrichment processes see ISchindler & Diaferid 
|2008| - Chapter 17, this volume). 

2 Sources of metals 

Most of the metals from O up to the Fe-group are produced by supernovae. The supemovae 
can be roughly divided into two groups: Type la supemovae (SN la) and core collapse su- 
pernovae (SNcc)- 



Table 1 Yields of elements relative to Fe in the Solar units of lGrevesse & Saliva 1 <1998h for different SN la 
models from Iwamoto et al. ( 1999). 



Element 


W7 


W70 


WDDl 


WDD2 


WDD3 


CDDl 


CDD2 





0.0314 


0.0280 


0.0214 


0.0136 


0.0105 


0.0235 


0.0114 


Ne 


0.00443 


0.00215 


0.00157 


0.000960 


0.00107 


0.00180 


0.000714 


Mg 


0.0221 


0.0390 


0.0219 


0.0109 


0.00584 


0.0231 


0.00976 


Si 


0.356 


0.317 


0.706 


0.453 


0.316 


0.746 


0.413 


S 


0.405 


0.410 


0.843 


0.544 


0.380 


0.890 


0.496 


Ar 


0.378 


0.433 


0.846 


0.555 


0.384 


0.896 


0.506 


Ca 


0.326 


0.460 


0.908 


0.604 


0.427 


0.972 


0.566 


Fe 


= 1.0 


= 1.0 


= 1.0 


= 1.0 


= 1.0 


= 1.0 


= 1.0 


Ni 


3.22 


2.18 


0.957 


1.25 


1.41 


0.914 


1.26 



SN la are most likely thermonuclear explosions of accreting white dwarfs. When the 
white dwarf reaches the Chandrasekhar limit, carbon ignition in the central region leads to a 
thermonuclear runaway. A flame front then propagates through the star at a subsonic speed 
as a deflagration wave. In the delayed-detonation models, the deflagration wave is assumed 
to be transformed into a detonation at a specific density. SN la produce a large amount of 
Fe, Ni, and Si-group elements (Si, S, Ar, and Ca). Contrary to SNcc, they produce only very 
small amounts of O, Ne, and Mg. In Fig.[T]and Table[T]we show the theoretically calculated 
yields for diff^erent SN la models (Iwamoto et al. 1999). We show the yields of elements 
relative to the yield of Fe in the Solar units of lOrevesse & Sauvall ([1998). Convective defla- 
gration models are represented by the W7 and W70 models. The WDDl, WDD2, WDD3, 
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Fig. 1 Yield s of elements relative to Fe in the Solar units of IGrevesse & Sauvall jl998l) for different SN la 
models from llwamoto et ai]U999l) . 



CDDl, and CDD2 models refer to delay ed-detonation explosion scenarios and the last digit 
indicates the deflagration to detonation density in units of 10^ g cm"^. The "C" and "W" 
refer to two different central densities (1.37 x 10'' and 2.12 x 10^ gem"-', respectively) in 
the model at the onset of the thermonuclear runaway. The relative yield of the Si-group ele- 
ments can be a good indicator of the incompleteness of the Si-group burning into Fe-group 
elements. Delayed detonation models have typically higher yields of Si-group elements rel- 
ative to Fe than the deflagration models. The Ni/Fe abundance ratio in the SN la ejecta is 
an indicator of neutron-rich isotope production, which depends on the electron capture ef- 
ficiency in the core of the exploding white dwarf The deflagration models produce larger 
Ni/Fe ratios than the delayed detonation models. 

In the left panel of Fig. [2] and Table [2] we show the theoretical yields for three SNcc 
models integrated over Salpeter and top-heavy initial mass fun ctions (IMFs) betw een 10 and 
50 M» llWooslev & Weave3ll995l : Iciiiefli & Limon gi 2004; Nomoto et al.ll2006h . See also 
iBorgani et al.ll2008l - Chapter 18, this volume for more details on stellar population models. 
In the right panel of Fig. |2] we show the SNcc yields for the Salp eter IMF for differen t 
progenitor metallicities (where Z = 0.02 refers to Solar metallicity; [ Nomoto et al] |2006l) . 
The Fe yield of the model by Woosley & Weaver ( 1995) is about an order o f magnitude 
lower than that in the models by Chieffi & Limongi (2004) and Nomoto et al.l 12006) . The 
plots also show that abundance ratios of 0/Mg and Ne/Mg can be used to discriminate 
between different IMFs. The SNcc yields for progenitors of different metallicities differ 
mainly in 0/Ne/Mg ratio. Because the progenitor metallicities and perhaps also the IMF 
evolve, the applicability of comparing the observed data with a single average set of yields 
may be limited. However, the abundances of these elements in the ICM may still be used 
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Fig. 2 Left panel: Yields of elements relative to Fe in the Solar units of iGrevesse & Sauva 1 (m^ for 
different SNcc models with a progenitor metallicity of Z = 0.004 integrated over Salpeter and top-heavy 
initial mass functions between 10 an d 50 Mp. Right panel : SNcc yields for different progenitor metallicities 
(Z = 0.02 refers to Solar metallicity; iNomoto et alj200"3) integrated over the Salpeter IMF. 



to infer to what extent was the dominant fraction of the SNcc progenitors preenriched, and 
what was its IMF. 

We note, that while the abundance s in Fig. [T] and [2] are shown with respect to the So- 
lar values of iGrevesse & Sauvall IT993) , in literature and also in several figures of this re- 
view the abundances are often shown with respect to the outdated Solar abundances by 
lAnders & Greves s? ("198^ or with respect to newer sets of Solar and proto-solar abun- 
dance s by Loddersi t2003i) . The more recent Solar abundan ce determinations of O, Ne , and 
Febv lLoddersl ( l2003l)" are ~30 % lower than those given by I Anders & Grevessel ( fl989l) .The 
Solar abundances of O and Ne reported by 'Grevesse & Sauval' (1998") are higher and the 
abundance of Fe is slightly lower than those reported by Lodders, C2003,) . A comparison of 
diff^erent Solar abundance values is shown in Table[3l 



Table 2 Yields of elements relative to Fe in the Solar units of'Grevesse & SauvaJ fT99i) for different SNcc 
models with a progenitor metallicity of Z = 0.004 integrated over Salpeter and top-heavy initial mass func- 
tions between 10 and 50 Mq. 



Element 


Nomoto et al. (2006) 
Salpeter top-heavy 


Chieffi & 
Salpeter 


Limongi (2004) 
top-heavy 


Woosley & Weaver (1995) 
Salpeter top-heavy 


O 


3.18 


10.8 


2.13 


6.67 


33.0 


79.6 


Ne 


3.55 


10.7 


2.53 


7.59 


15.4 


31.8 


Mg 


2.71 


6.63 


1.80 


3.94 


13.8 


28.5 


Si 


2.88 


6.50 


1.82 


3.55 


29.9 


50.9 


S 


1.95 


4.49 


1.33 


2.72 


23.9 


42.7 


Ar 


1.87 


4.24 


1.49 


3.08 


30.5 


57.1 


Ca 


1.63 


3.70 


1.37 


2.80 


22.2 


31.1 


Fe 


=1.0 


=1.0 


= 1.0 


= 1.0 


= 1.0 


= 1.0 


Ni 


1.10 


0.475 


1.11 


0.501 


497. 


208. 



While elements from O up to the Fe-group are produced by supemovae, the main sources 
of carbon and nitrogen are still being debated. Both elements are believed to originate from 
a wide range of sources including winds of short-lived massive metal rich stars, longer- 



5 



Table 3 Comparison of Solar abundances by Anders & Grevesse (AG; Il989h . iGrevesse & Sauva] (GS; 
119981) . and of the prot-solar abundances bV iLodders, C20Q3.) on a logarithmic scale with H=12. These abun- 
dances are used at several places of this volume. 



element 


AG 


GS 


Lodders 


H 


12.00 


12.00 


12.00 


He 


10.99 


10.93 


10.98 


C 


8.56 


8.52 


8.46 


N 


8.05 


7.92 


7.90 





8.93 


8.83 


8.76 


Ne 


8.09 


8.08 


7.95 


Mg 


7.58 


7.58 


7.62 


Si 


7.55 


7.55 


7.61 


S 


7.21 


7.33 


7.26 


Ar 


6.56 


6.40 


6.62 


Ca 


6.36 


6.36 


6.41 


Fe 


7.67 


7.50 


7.54 


Ni 


6.25 


6.25 


6.29 



lived low- and in tennediate-mass stars, and also an early generation of massive stars (e.g. 
iGustafsson eTal. 1999; Chiappini et al. 2003; Meynet & Maeder 2002). 

In the Galaxv. lShi et alJ ll2002li found that C is enriched by winds of metal-rich massive 
stars at the beginning of t he Galactic disk evo l ution, while at a later stage it is produced 
mainly by low-mass stars. iBensbv & Feltzind ( l2006h found that the C enrichment in the 
Galaxy is happening on a time scale very similar to that of Fe. They conclude that while in 
the early Universe the main C contributors are massive stars, C is later produced mainly by 
asymptotic giant branch (AGB) stars. 

Nitrogen is produced during hydrogen burning via the CNO and CN cycles as both a 
primary and secondary element. In primary nucleosynthesis its production is independent 
of the initial metallicity of the star. Primary produc tion of N happens during hydrogen shell 
burning in intermediate mass stars of ~ 4 - 8 Mq jMatteucci & Tosilll985h Stellar models 
that include the effe cts of rotation indicate th at massive stars between 9 and 20 Mq may also 
produce primary N ( iMaeder & Mevnetll2000h . In secondary production, which is common to 
stars of all masses, N is synthesised from C and O and its abundance is therefore proportional 
to the initial metallicity of the star. The primary versus secondary production of N can be 
studied by investigating the N/0 ratio as a function of 0/H ratio. In the case of primary N 
production the N/0 ratio will be constant. For the secondary production we will observe 
a linear correlation between the logarithms of N/0 and 0/H. The combination of primary 
and secondary nucleosynthesis will produce a non-linear relation. The ICM abundance of 
C and N was measured only in a few bri" ht nearby galaxy clusters and elliptical galaxies 
jPeterson et alj[2003l ; IWerner et al.ll2006a ). However, future instruments (see Sect. [8ll will 
allow to measure the abundances of these elements in many clusters. This will be important 
for a better understanding of the stellar nucleosynthesis. 



3 Abundance studies before XMM-Newton and Chandra 



The discovery of the Fe-K line emis sion in the spectrum o f the Perseus cluster by the Arie7 V 
satellite (see the left panel of Fig. [3l lMitchell et"ai]|l976h and in Coma and Virgo by OSO-8 
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Fig. 3 Left panel: Spectrum of the Perseus cluster obtained by Ariel 5. The first cluster spectrum with ob- 
served Fe-K line emission (hump at ~7 keV). From Mitchell et al. ( 1976). Right panel: Line em ission ob- 
served in M 87 with the Einstein satellite. From iSarazin l ,1988i\ based on the work published bv lLea et alj 
fT982h. 



dSerlemitsos et al.lll977h confirmed that the X-ray emission of galaxy clusters is predomi- 
nantly thermal radiation from hot intra-cluster gas rather than inverse Compton radiation. 
These observations showed that the hot plasma in clusters of galaxies contains a signifi- 
cant portion of processed gas, which was ejected from stars in the cluster galaxies. Sub- 
sequent spectroscopic analysis of cluster samples observed with OSO-8 and HEAO-1 A2 
satellites revealed that the ICM has an Fe abundance of about one-third to one-half of the 
Solar value (Mushotzky et al. 1978; Mushotzky 1984). Combining spectra obtained by Ein- 
stein and Ginga. iWhite et all | |1994|) found an indication of a centrally enhanced metallicity 
in four cooling flow clusters. Unfortunately, the spectrometers on the Einstein observatory 
were not sensitive to the Fe-K line emission, because the mirror on the satellite was not sen- 
sitive to photon energies above ~4 keV. However, the Solid-State Spectrom eter (SSS) and 
the h igh-resolution Focal Plane Crystal Spectrometer (FPCS) on Einstein (Giacconi et al.l 
Il979h allowed to detect emission lines at low energies. Using the SSS the K lines from Mg, 
Si, and S and the L lines from Fe were detected in the spectrum of M 87 (see the ri ght panel 
of Fig.©, Perseus A 49 6, and A 576 dUea et al. 1982: Mushotzk y et a l. 1981: N ulsen et all 
[l982; Mus hotzky! 1 19841) . Using the FPCS spectra of M 87 the O vm Kor line was detected 
ICanizares et al.lll979h . implying an 0/Fe ratio of 3-5, and the relative strength of various 
Fe-L line blends showed that the gas cannot be at a single temperature. 

However, until the launch of ASCA in 1993, Fe was the only element for which the 
abundance was accurately measured in a large number of clusters. ASCA allowed to detect 
the emission features from O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni in the spectra of a number of 
clusters. Furthermore, it allowed to accurately determine the Fe abundances out to redshift 
z » 0.5. ASCA data revealed a lack of evolution in the Fe abundance out to redshift z ~ 
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Fig. 4 Left panel: The radial distribution of the logarithm of the Si/Fe ratio expressed in Solar units (scale 
of Anders & Grevesse 1989) for a sample of clusters observed with ASCA (from Finoauenov et al. 2 OO0I). 
Right panel: The average radial distribution of the Fe abundance in Solar units (scale of Grevesse & Sauval 
[1998) for a sampl e of cooling core (fille d circles) and non-cooling core (empty circles) clusters observed wfith 
BeppoSAX (from lDe Grandi^et al.ll2004ft . 



0.4 ('Mushotzky & Loewenstein"l997'; 'Rizza e t aDll998h and no evidence for a decrease at 
higher redshifts (Donahue et al. 1999). ASCA data we re the first to show th at in cooling 
core clusters the metallicity increases toward the ce ntre 1 Fukazawa et al.lll994) . 

By analysing data of four clusters of galaxies, iMushotzkv et al.l ( Il996h found that the 
relative abundances of O, Ne, Si, S, and Fe are consistent with an origin in SNcc- They 
suggest that SNcc provided a significant fraction of metals in the ICM. Ipukazawa et al.l 
mm showed that the Si abundance and the Si/Fe ratio increases from the poorer to the 
richer clusters, suggesting that the relative contribution of SNcc increases toward richer 
clusters. They prop ose the possibility t hat a considerable fraction of SNcc products escaped 
the poorer systems. [Finoguenov et alj ( I200Q) used ASCA data to show that SN la products 
dominate in the cluster centres and the SNcc products are more evenly distributed (see the 
left panel of Fig.|4]l. 

However, the large and energy dependent point-spread function (~2' half-power ra- 
dius) of ASCA did not allow to investigate the spatial abundance distributions in detail. 
iDe Grandi &Molendil l l200 ll) took advantage of the better spatial resolution (~r half-power 
radius) of BeppoSAX and they measured the radial Fe abundance profiles for a sample of 17 
rich nearby clusters of galaxies. They found that while the eight non-cooling core clusters 
in their sample have flat Fe abundance profiles, the Fe abundance is enhanced in the centra l 
regions of the cooling core clusters (see the right panel of Fig. l4ll. loe Grandi et al.l ( |2004|) 
show that the mass associated with the abundance excess found at the centre of cool core 
clusters can be entirely produced by the brightest cluster galaxy (BCG). 

The bandpass of the i?OSAT PSPC was well suited for X-ray bright galaxies and groups 
of galaxies with temperatures of ~ 1 keV and its PSF allowed for spatially resolved spectral 
analysis on a half-arcminute scale. Buote (2000a) investigated the deprojected abundance 
gradients in 10 galaxies and groups. In 9 of the 10 systems they found an abundance gradient, 
with the Fe abundance one to several times the Solar value within ~10 kpc, decreasing to 
-0.5 Solar at 50-100 kpc distance. 

The emerging picture from the abundance studies with ASCA and BeppoSAX data was 
that of an early homogeneous enrichment by SNcc> which produce c-elements in the pro- 
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Fig. 5 These plots illustrate the influence of the temperature model on the best-fit abundances. Left panel: 
Radial Fe a bundance distribution in NGC 1399 modelled with single- and two-temperature models (from 
lBuotel20"o3) . Right panel: Radial S abundance distribution in the cluster 2A 0335+096 modelled with single , 
two, and multi-temperature models (triangles, squa res, and circles, respective ly; from lWemer et al.ll2006bl) . 
Abundances in both figures are in the Solar units of lAnders & Grevessd jl989h . 



tocluster phase, and a subsequent more centrally peaked enrichment by SN la which have 
longer delay times and continue to explode in the cD galaxy for a long time after the cluster 
is formed. 

Analysing an d stacking all the d a ta of clusters of galaxies observed with ASCA into 
temperature bins, *Baumgartner et al. 1 12005") found a trend showing that the relative con- 
tribution of SNcc is larger in higher temperature clusters than in the clusters with lower 
temperature. Moreover, they found that Si and S do not track each other as a function of 
temperature, and Ca and Ar have much lower abundances than expected. They conclude 
that the ff-elements do not behave homogeneously as a single group and the trends indicate 
that different enrichment mechanisms are important in clusters of different masses. 



4 The possibilities and limitations of elemental abundance determinations 

4. 1 Spectral modelling and the most common biases 

The correct modelling of the temperature structure is crucial for the elemental-abundance 
determinations. Fitting the spectrum of a multi-temperature plasma with a simple single- 
temperature model results in best-fitting elemental-abundances that are too low. This effect 
is the most important in galaxies, groups of galaxies and in the cooling cores of clusters. 
Many studies have found significantly subsolar values for the Fe abundance in groups of 
galaxies, with abundances generally less than the stellar values in these systems, and less 
than those observed in galaxy clusters. To reconcile these observations with the chemical en- 
richment models, it has been postulated that galaxy grou ps accrete primordial gas after they 
have been expelling gas during most of their evolution jRenzinilll997h . iBuote & CanizaresI 
fl994) showed that if a galaxy spectrum has intrinsically two temperature c omponents, the 
single-temperature model can give a significantly underestimated metallicitv. iBuotd ( l2000bl) 
demonstrated that in a spectrum characterised by two temperatures, one below 1 keV and 
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one above 1 keV, with an average of ~ 1 keV, the lower temperature component preferen- 
tially excites emission lines in the Fe-L complex below 1 keV (Fe xvii-xxi), while the higher 
temperature component excites the Fe-L lines from ~1-1.4 keV (Fexxi-xxiv). If these com- 
ponents contribute approximately equally to the emission measure, the spectral shape of the 
Fe-L complex will be broader and flatter than the narrower Fe-L complex produced by a 
single-temperature model with the average temperature of the two components. In order to 
fit the flat spectral shape at 1 keV, the single-temperature model will fit a lower Fe abundance 
and will increase the thermal continuum. The strongest bias is at the low temperatures, for 
higher temperature plasmas the bias is smaller as the Fe-K line-e mission becomes more im- 
portant for the determination of the Fe abundance. iBuotd (l2000b l) also shows an excess in the 
Si/Fe and S/Fe ratios inferred for two-temperature plasmas fit by single-temperature mod- 
els. He notes that the lower temperature component in the two-temperature plasma produces 
stronger Si and S emission, and the single-temperature model underestimates the contin- 
uum at higher energies. Therefore the Si/Fe and S/Fe ratios of multi-temperature IC M in- 
ferred u sing single-temperature models will yield values in excess of the true values. ISuotd 
shows that this effect is the most important for the hotter clusters of galaxies. This 
bias may have played a role in the excess Si/F e ratios and trends derived using ASCA data 
jMushotzkv et al.lll99d iFukazawa et al.lll998l ; iFinoguenov et al.ll2000l ; iBaumgartner et al.l 
l2005b . 

Several authors (e.g. iBuote et al.ll2003l : [Wemer et al]|2006bl : [Matsushita et aDl2007ah . 
using XMM-Newton and Chandra data, showed that the best-fitting abundances in the multi- 
phase cooling core regions of clusters and in groups increase when two-temper ature models 
or mo re complicated difi^erential emission measure (DEM) models are used. iBuote et aP 
( l2003l ) show that the abundances of O, Mg, Si, and Fe in the centre of NGC 5044 obtained 
using a two-temperature model are higher than those obtained using a single-temperature 
model. They also show that the Si/Fe ratio in the centre is significantly lower when a two- 
temperature model is used. Similar results were obtained for NGC 1399 (see the left panel 
of Fig. [5] Buote 200 2). Analysing deep XMM-Newton data of the cluster 2A 0335-H096, 
IWemer etall J2006bl ) showed that the abundances of all fitted elements (Mg, Si, S, Ar, Ca, 
Fe, Ni) in their model increased as they went from single- to two-temperature model, and 
to the multi-temperature DEM model. However, not only their absolute values increased 
but also the Si/Fe and the S/Fe ratios changed. iMatsushita et al] ( 12007 af ) show an increase 
from single- to two- to multi-temperature models for O, Mg, Si, and Fe in the core of the 
Centaurus clus ter. In the right pan el of Fig. [5] we show the radial S profile obtained for 
2A 0335H-096 jWemer et al.ll2006bh with single-, two-, and multi-temperature models. 

As discussed in i Buotd ( l2000al) and Buot e et al. I ( l2003h in the lower mass systems (~ 
1 keV) the best fit Fe abundance is sensitive to the lower energy limit used in spectral fitting. 
If the continuum at the low energies is not properly modelled, the equivalent widths of the 
lines in the Fe-L complex (and the o xygen line em ission) will be inaccurate and the abun- 
dance determination will be incorrect. iBuote et al.l (12003) show the best-fitting Fe abundance 
in NGC 5044 determined with the lower energy limit at 0.7, 0.5, and 0.3 keV. They show 
that the results for each energy limit are significantly diff^erent and they obtain the highest 
Fe abundance for the lowest energy limit, which allows them to properly fit the continuum. 

Since the ~ 0.2 keV component of the Galactic foreground emits O viii line emission, 
the uncertainties in the foreground/background model can cause systematic uncertainty in 
determining the oxygen abundance in the ICM. The temperature and flux of the ~ 0.2 keV 
Galactic foreground component is position dependent on the sky and can not be properly 
subtracted with blank fields (for detailed discussion see lde Plaa et al.ll2006l) . 
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Systematic uncertainties in the abundance determination can be introduced if bright 
point sources are not excluded from the spectral fits, especially in the lower surface bright- 
ness regions of clusters. The contribution of point sources raises the continuum level, sys- 
tematically reducing the observed line equivalent widths. If their contribution is not properly 
modeled, the best fitting thermal model will have a higher temperature and lower elemental 
abundances. 



In general tw o plasma mod els are used for fitting cluster spectra: MEKAL I Mewe et al 



I9851 . 1 19861 : Ixaastra 1992 : Liedahl et al.ll 19951 : iMewe et al.|[r995) and APEC jSmith et al 



20011) . These codes have differences in their line libraries and in the way of modelling the 



atomic physics. This allows to test the robustness of the abundance determinations. System- 
atic effects on the best fit abund a nces due to the adopted pl asma model were investiga ted for 
example bv lSuoteet alj ( l2003h , ISanders & FabianI ( l2006h . and Ide Plaa et all i200% . They 
conclude that despite small differences in the determined temperature structure, the derived 

abundances are generally consistent v yithin the errors . 

Using simulated galaxy clusters iKapferer et all ( l2007ah show that the more inhomo- 
geneously the metals are distributed within the cluster, the less accurate is the observed 
metallicity as a measure for the true metal mass. They show that in general the true metal 
mass is underestimated by observations. If the distribution of metals in the ICM is inhomo- 
geneous, the true metal mass in the central parts (r < 500 kpc) of galaxy clusters can be up 
to three times higher than the metal mass obtained by X-ray observations. The discrepancies 
are due to averaging. As the metallicity and temperature are not constant throughout the 
extraction area, the integration of thermal Bremsstrahlung and of line radiation can lead to 
underestimated metal masses. 



4.2 Possibilities and limitations of current instruments 

XMM-Newton with its large effective area and superb spectral resolution, Chandra with 
its unprecedented spatial resolution, and Suzaku with its low background and good spec- 
tral resolution at low energies largely enhanced our knowledge on chemical abundances in 
galaxies, groups, and clusters of galaxies. In this section we describe the possibilities and 
the systematic uncertainties of abundance studies with these instruments. 

4.2.1 XMM-Newton 

Due to its large eff"ective area and good spectral resolution, XMM-Newton is currently the 
best suited instrument to study abundances in the ICM. XMM-Newton has three identical 
telescopes and two sets of detectors on board. The first set of instruments is the European 
Photon Imaging Camera (EPIC), consisting of two metal oxide semiconductor - MOS 1 and 
M0S2 - CCD arrays collecting 50 % of light from two telescopes, and the EPIC-pn CCD ar- 
ray collecting 100 % of light from one telescope on board. The second set of instruments are 
the Reflection Grating Spectrometers (RGS) that collect 50 % of light from two telescopes. 
Both sets of instruments play an important role in the abundance studies. 

The advantages of EPIC are its broad bandpass, and good spatial and spectral resolu- 
tions, that allow us to detect the emission lines from O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni 
(see the top panel of Fig. [S}. Furthermore, for deep observations of bright clusters, EPIC 
allows us to derive the spatial distribution of elements (see Sect.]?)- 

Unfortunately, the abundances determined in the cluster outskirts with low surface bright- 
ness using XMM-Newton remain uncertain due to the high background level. The value 
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Fig. 6 Top panel: The line spectrum of the cluster 2 A 0335+096, as observed with XMM-Newton EPIC 
(fromfWemer et al. 2006b). Bottom panel: Line spectrum of M 87. as observed with XMM-JVewton RGS 
(from Werner et al. 2006a). After the fitting, the line emission was set to zero in the multi-temperature plasma 
model indicated on the y-axis. While the CCD spectra are plotted as a function of the observed energy, the 
grating spectra are shown as a function of the observed wavelength. 
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of the best fit te mperature in these r egions is very sensitive to the level of the subtracted 
background (e.g. Ide Plaa et al.ll2006[) and as we show above in Sect. 14. l[ the temperature 
modelling is very important for the correct abundance determination. 

The O line emission in the spectrum is strong, but its abundance is difficult to measure 
because it is sensitive to the assumed model for the Galactic foreground (see Sect. 14.1b 
and there are still remaining calibration uncertainties at low energies. The ls-2p Ne lines 
at 1.02 keV are in the middle of the Fe-L complex (lying between 0.8 to 1.4 keV). The 
resolution of the EPIC cameras is not sufficient to resolve the individual lines in the Fe-L 
complex, which makes the Ne abundance determination by EPIC uncertain. However, the 
abundances of O and Ne in the cores of cooling core clusters, observed with XMM-Newton 
with sufficient statistics, can be well determined using the high spectral resolution of the 
RGS. 

The systematic differences in the elfective-area between EPIC-MOS and EPIC-pn are 
of the order of 5-10 % in certain bands ( Kirsch 2006). By using data with high photon 
statist ics and fitting spectra obtained by EPIC-MOS and EPIC-pn separately, Ide Plaa et alj 
l l2007h investigated the magnitude of the systematic errors on the best-fitting abundances 
of elements heavier than Ne due to the uncertainties in the effective area calibration. They 
found that the elemental abundances most inffuenced by the calibration uncertainties are Mg, 
Si, and Ni. de Plaa et al., (2007) show that the systematic uncertainties on Mg are the largest, 
since the eff^ective-area of EPIC-pn has a dip with respect to the eff^ective-area of EPIC-MOS 
at the energies of the magnesium lines. This causes the Mg abundance derived by EPIC-pn 
to be underestimated. For the abundances of Si and Ni they estimate that the systematic 
uncertainties are 1 1 % and 19 %, respectively. Although EPIC-MOS and EPIC-pn show 
large systematic differences at the energy of the Mg lines, the Mg abund ances derived by 
EPIC-MOS an d Chandra ACIS for the Centaurus cluster are consistent ( iMatsushita et al.l 
l2007d : ISanders & Fabiaj|200d) , and the Mg abundances derived by the EPIC-MOS and 
RGS using the deep observations of 2A 0335H-096 also agre e with each other. 

The S abundance determination with EPIC is robust (e.g. Ide Plaa et al .120071) . Ar and Ca 
abundances should also be accurately determined, however their signal is usually weak, so 
the level of their observed equivalent widths may be somewhat sensitive to the continuum 
level. 

The RGS has a high spectral resolution but only a limited spatial resolution in the cross- 
dispersion direction. It allows for relatively accurate abundance measurements of O, Ne, 
Mg, Fe, and in the case of deep observations of nearby bright clusters and elliptical galaxies 
even the spectral lines of C and N are detected (see the bottom panel of Fig.|6ll. In the case of 
deep observations of nearby bright clusters, RGS allows us to extract spectra from typically 
five bins that are ~r wide in the cross-dispersion direction and the photons in each bin are 
collected from ~10' long regions in the dispersion direction. This way RGS allows us to 
determine projected abundance profiles. 

The spectral resolution of the RGS is, however, limited by the broadening of the ob- 
served line profiles caused by the spatial extent of the source along the dispersion direction. 
This makes RGS spectroscopy useful only for cooling core clusters with strongly peaked 
surface brightness distribution. In order to account for the line broadening in the spectral 
modelling, the model is convolved with the surface brightness profile of the source along 
the dispersion direction. Because the radial profile for a spectral line produced by an ion can 
be different from the radial surface brightness profile in a broad band, in practice the line 
profile is multiplied with a scale factor, which is a free parameter in the spectral fit. Different 
best fit line profiles for the Fe li nes and for the O li ne emission can suggest a different radial 
distribution for these elements Ide Plaa et al.ll200^ . 
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4.2.2 Chandra 

Although the Advanced CCD Imaging Spectrometer (ACIS) on Chandra has a much smaller 
effective area than the EPIC on XMM-Newton, which makes it less suitable for abundance 
studies of the ICM, it is the instrument of choice for studies of metal enrichment in galaxies, 
where a good spatial resolution is required. Chandra is able to resolve the substantial fraction 
of X-ray binaries in galaxies as point sources, so their contribution can be excluded from the 
spectral modelling of the inter-galactic medium (IGM). The excellent spatial resolution also 
makes it possible to investigate the spatial temperature variations, which can bring a bias into 
the determination of abundances (see Sect. l4.lt . The spectral resolution of ACIS is sufficient 
to reliably determine the abundances of a number of elements in the ICM and IGM. With 
deep enough exposure times Chandra ACIS can be iised to study the ICM abiindanc es of O, 
Ne, Mg, Si, S, Ar, Ca, Fe, and Ni (e.g. lSanders et al.ll2004l ; [Sanders & Fab"iar]|2006l) . 

iHumphrev et al.l ( l2004h discuss the systematic uncertainties in the abundance determi- 
nations due to the uncertainties in the absolute calibratio n of ACIS. After corre cting for the 
charge transfer inefficiency (CTI) using an algorithm bv lTownslev et al.l ( l2002h . they found 
a change of zlZpe - 0.3 in the best fit Fe abundance of NGC 1332. 

I Sanders & FabianI ( l200d) investigate the systematic uncertainties on the abundance de- 
terminations with the S3 chip of Chandra ACIS in the Centaums cluster. They conclude that 
for the best determined elements the non-CTI-corrected Chandra and XMM-Newton results 
agree well. The only element where the Chandra and XMM-Newton data do not agree is 
oxygen. There has been a buildup of contamination on the ACIS detectors, which influences 
the effective area calibration at the low energies making the O abundance determination dif- 
ficult. The iridium coating of the Chandra mirrors produces a steep drop in the effective area 
at around 2.05 keV, introducing a potential systematic uncertainty in the Si abundance de- 
termi nation. However, the Si abundances determined by Chandra and XMM-Newton agree 
well jSanders & Fabiaj|2006h . 



4.2.3 Suzaku 

The advant age of the X-ray Im aging Spectrometer (XIS) detectors jKovama et alj|2007l) 
on Suzaku jMitsuda et al.ll2007b compared to EPIC on XMM-Newton is its better spectral 
resolution below 1 keV, and the low background, providing a better sensitivity at the O 
line energy. Suzaku also allows for better measurements of the Mg lines in the low surface 
brightness cluster outskirts, which is difficult with EPIC/MOS due to the presence of a strong 
instrumental Al line close to the expected energy of the observed Mg lines. 

The relatively broad point spread function of XIS with a half-power diameter of ~2' 
jSerlemitsos et al.ll2007h does not allow detailed investigations of the spatial distribution of 
elements. 

The effective area of the XIS detectors below 1 keV is affected by carbon and oxygen 
contamination of the optical blocking filter. The contamination is time and position de pen- 
dent. It s column density and spatial distribution are still not well known. Matsushita e t al.l 
l l2007bl) estimate that the uncertainty in the thickness and chemical composition of the con- 
taminating column results in a ~20 % systematic error on the O abundance. 
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Fig. 7 Left panel: Radial distribution of Fe and O in M 87 determined using XMM-Newton RGS (from 
IWemer e t al. 2006a). The abundances are shown relative to the proto-solar values given by Lodders i2C^. 
Right panel: Radial distribution of the O/Fe ratio determined for a sample of clusters using XMM-Newton 
EPIC (from jT amura et al. 2004J. The abundance ratios are shown relative to the Solar values given by 
lAnders & Grevessel 11989!) . 



5 Spatial distribution of elements 

5.1 Radial abundance profiles 

Observations with XMM-Newton and Chandra confirmed the centrally peaked metallic- 
ity distribution in the cooling core clusters and the flat distri bution of metals in the non- 
cooling core clusters dVikhlinin et"d]|2005l : IPratt et al.ll2007h previously seen with Bep- 
poSAX. XMM-Newton and Chandra also confirmed the strong contribution of SN la to 
the enrichment of cluster cores found by ASCA. The observations are consistent with the 
picture of an early contribution by SNcc and later enrichment of cluster cores by SN la. One 
of the first abundance studies with XMM-Newton, performed on the cooling core cluster 
Abell 496, showed a radially flat distribution of O, Ne, and Mg and centrally peaked distri- 
bution of Ar, Ca, Fe, and Ni jXamura et al The abundance peak in the cluster core is 
consistent with the idea that the excess metals were produced by SNe la in the cD galaxy. 
However, the strongly increasing Si/Fe ratio toward the outskirts seen with ASCA was not 

confirmed with XMM-Newton. 

An alysing XMM-Newton EPIC data of M 87. lB6hringer et aD(l200lh and lFinoguenov et al.l 
( l2002h found rather low O/Fe rati o in the cluster core, that incre ased toward the outer re- 
gions. Analysing the same dataset, iGastaldello & Molendil ( l2002h questioned these results, 
showing that the O/Fe ratio in th e inner 9' is c o nstant. As shown in the left panel of Fig. |7] 
usme XMM-Newton RGS data, IWemer et al] ( l2006ah found a centrally pe aked Fe abun- 
dance g radient and a flat O distr ibutio n in M 87 confirming th e results of Bohringer et alj 
bOOlh . iFinoguenov et al.l (l2002h . and [Matsushita et alj ( l2003h (the O abundance profiles 



found bv' Gastaldello & Molendl ilOoj) were dififerent because of calibration problems at 
low energies in the early days of the mission). The RGS data reveal an O/Fe ratio of 0.6 
in the cluster core and a ratio of 0.9 at a distance of 3' (14 kpc) from the centre. However, 
the XMM-Newton EPIC data show that the Si abundance has a similar gradient in M 87 
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as the Fe abundance. A similar abundance pattern with a rather flat radial distribution of O, 
and centrall y peaked Si, and Fe abundances is observed in the cluste r sample analysed b y 
IXamura eTa l. ( 2004, s ee the right panel of Fig.|7]i, in t he Perseus cluster I Sanders et al.l2004[) . 



in Abell 85 (Durret eTaTlbOOSt) , in Sersic 159-03 ( de Plaa et al 



(^^ mer et al. ,2006h') . in the Cent aums cluster (Matsushita et al 



20061} in 2A 0335+096 



2007al : ISanders & FabianI 



l2006h . and in the Fornax cluster ( iMatsushita et alJ 2007a[ ). The abundances of Si and Fe 



have similar gradients also in the group of galaxies NGC 5044 (' Buote et al.ll2003h , with a 
constant radial distribution of Si/Fe= 0.83. The radial distribution of O in this group is con- 
sistent w ith being flat. Suz aku observations of the clusters A 1060 dSato et al.ll2007ch and 
AWM 7 dSato et alJl2007ah also showed rather low O abundance of around 0.5 Solar within 
about 5' from the core without a clear abundance decline in the outer regions. 

The radial metallicity gradients in cooling flow clusters often display an inversion in 
the centre (e.g. ISanders & Fabianll2002l) . This apparent metallicity drop in the very central 
region is often the result of an oversimplified model in the spectral analysis (see Sect.|4]and 
iMolendi & GastaldelldlioOlh . However, in some cases the metallicity drop in the cluster 
core is robust. In the Perseus cluster the abundances peak at the radius of 40 kpc and they 
decrease at smaller radii. This inversion does not disappear when extra temperature compo- 
nents and po wer-law models are included, or when projection eff'ects are taken into account 
dSanders & Fabian,2007). 



iBohringer et al. found that long enrichment times (>5 Gyr) are necessary to pro- 



duce the observed central abundance peaks. Since the classical cooling flows, or a strongly 
convective intra-cluster medium, or a complete metal mixing by cluster mergers would de- 
stroy the observed abundance gradients, they conclud e th at cooling cores in cluste rs are 
preserved over very long times. lBShringer et al.l ll2004h and [Matsushita et aP ( l2007bh show 
that the iron-mass-to-light ratio in the central region of the Centaums cluster (the mass of 
Fe relative to the B-band luminosity of the stellar population of the cD galaxy) is much 
higher than that of M 87. This indicates that the accumulation time scale of the SN la prod- 
ucts was much higher in Centaurus than in M 87. In agreement with the conclusions of 
iDe Grandi et al. I (l2004l) they also conclude that the innermost part of the ICM is dominated 
by gas originating mainly from the stellar mass loss of the cD galaxy. 

If the metals originate from the stellar population of the cD galaxy, then in the absence 
of mixing the abundance profiles would follow the light profiles. The observed metal profiles 
are much less steep than the lig ht profiles, which suggests th at the injected metals are mixed 
and they diffuse to larger radii. iRebusco etal.lll2005l . l200^ model the diffusion of metals 
by stochastic gas motions after they were ejected by the brightest galaxy and compare their 
results with the peaked Fe abundance profiles of 8 groups and clusters. They suggest that 
the AGN/ICM interaction makes an important contribution to the gas motion in the cluster 
cores. 

In order to determine the metallicities of the stellar populations of galaxies and star 
clusters it has become a common practice to use the feature made up of the Mg b line and 
the Mg H band around 5200 A, which is known as the Mg2 index. The Mg2 index appears 
to be a relatively good metallicity indicator. In M 87, Centaurus, and the Fornax clusters 
the O and Mg abundances in the centre a re consisten t with the optical metallicity of the cD 
galaxy derived from the Mg2 index IMatsushita et al.t200 3. 2007b, aj. The Mg2 index mainly 
depends on the Mg abundance, but also on the total metallicity where the O contribution 
matters most, and depends weakly on the age population of stars. The agreement between 
the O and Mg abundances in the ICM with the stellar metallicities supports the view that the 
ICM in the cluster centres is dominated by the gas ejected from the cD galaxies. However, 
the Mg/Fe ratio in the stellar populations of elliptical galaxies is significantly higher than 
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that in the hot sas (iHumphrev & Buotell2006h . The observed difference shows that since the 
time when the current stellar populations of these ellipticals formed, a significant number of 
SNe la enriched the gas in addition to the stellar mass loss. 

Observations indicate that at large radii clusters have a flat abundance distribution of 
about 0.2 Solar. By using Suzaku data, iFuiita et all ( l2007h found that the gas is enriched to 
a metallicity of ~0.2 Solar out to the virial radii of clusters A 399 and A 401. These clusters 
are at the initial stage of a merger and their virial radii partially overlap. The authors interpret 
their finding as a result of early enrichme nt by galactic supe rwinds in the proto-cluster stage. 
However, the extraction region used by iFuiita et aD l l2007h is large and their measured Fe 
abundance might be significantly biased toward the value at about half of the virial radii of 
the clusters. 



5.2 2D distribution of metals 



Deep observations of bright clusters of galaxies with XMM-Newton and Qiandra allow 
to map the 2D distribution of metals in the ICM jSanders et alj|2004 iDurret et al.ll2005 ; 



O' Sullivan et al. 2005; Sauvageot et al. 2005; Finoguenov et al.ll2006l ; IWemer et al.ll2006h 



Sanders & Fabian 2006 ; Bagchi et al. .2006.) . In general the distribution of metals in clusters 



is not spherically synnmetric, and in several cases it shows edges in the abundance distribu- 
tion or several maxima and complex metal patterns. The 2D distribution of metals is a good 
diagnostic tool, which can be used together with the maps of thermodynamic properties to 
investigate the comp lex structure of cluster cores and the merging history of clusters (e.g. 
iKapferer et ai]|2006h . 

M 87 provides us with the unique opportunity to study in detail the role of the AGN 
feedback in transporting and distributing the metals into the ICM. The high X-ray sur- 
face brightness and small distance of M 87 allow a detailed study of the 2D distribution 
of several elements (O, Si, S, and Fe). By analysing deep (120 ks) XMM-Newton data, 
ISimionescu et al found that the spatial distribution of metals shows a clear enhance- 

ment along the radio lobes, where the rising radio emitting plasma ejected by the AGN 
uplifts cooler gas from the core of the cluster (see Fig.[8ll. Furthermore, they showed that the 
metallicity of the cool uplifted (kJ < 1.5 keV) gas is ~2.2 times Solar and the abundance 
ratios in and outside the multi-phase region associated with the uplifted gas are very similar. 
They estimate that the mass of the cool gas is »5 x 10^ Mq and it probably originates from 
the stellar winds enriched with SN la products. Approximately 2 x 10** yr are required to 
produce the observed mass of cool gas, indicating that the uplift of cool gas by AGN radio 
bubbles is a rare event. 



6 Element ratios and their reconstruction with supernova models 

The chemical elements that we see in the ICM are the integral yield of aU the diff^erent super- 
nova types that have left their specific abundance patterns in the gas prior and during cluster 
evolution. Since the launch of ASCA, which for the first time allowed the determination of 
abundances of elements other than Fe in the ICM (see Sect. [3]l, several groups tried to use 
the ICM abundance patterns to constrain the contribution of different kinds of supernovae to 
the ICM enrichment and to put constraints on the theoretical supernova models. 

In the following subsection we review these efforts. The still somewhat ambiguous re- 
sults of these efforts are summarised in section 6.2 
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Fig. 8 Map of the Fe abundance in M 87. The colour bar indicates the Fe abundance with respect to the 
proto-solar value of Lodders (2003), with red showing higher and blue lower abundances. Contours of the 
90 cm radio emission ( Owen et al...2000.) are overplotted. The blacli circle indicates the half light radius of 
M 87. Beyond the expected radial grad ient, one clearly sees the enhanced Fe abundance in the radio arms, 
especially within the Eastern arm. From lSimionescu et ai] )2007l) . 



6. 1 Results on reconstructed supernova models 

iMushotzkv eTaP ( Il996h investigated the elemental abundances in four clusters, removing 
the central regions, and concluded that the abundance ratios are consistent with the contri- 
bution of SNcc- Dupke & White (2000) found that most of the Fe (~ 60 - 70 %) in the 
cores of their investigated clusters comes from SN la, and they used these central regions 
dominated by SN la ejecta to discriminate between competing theoretical models for SN la. 
They showed that the observed high Ni/Fe ratio of ~4 is more consistent with the W7 con- 
vective deflagration SN la model than with the delayed detonation models. Stacking all the 
cluster data in the ASCA archive into temperature bins, Baumgartner et al ] (f2005) showed 
that the abundances of the most well determined elements Fe, Si, and S do not give a con- 
sistent solution for the fraction of material produced by SN la and SNcc at any given cluster 
temperature. They concluded that the pattern of elemental abundances requires an additional 
source of metals other than SN la and SNcc, which could be the early generation of Popu- 
lation III stars. 

The ASCA estimates were largely based on the Si/Fe ratio, which has the disadvantage 
of Si being contributed by both SN la and SNcc- The much improved sensitivity of XMM- 
Newton to measure the spectral lines of more elements at a higher statistical significance 
off ers be tter opportuni ties to constrain the supernova models. 

iFinoguenov et al.i ( I2OO2I) used XMM-Newton EPIC to accurately measure the abun- 
dances of O, Ne, Mg, Si, S, Ar, Ca, Fe, and Ni in the hot gas in two radial bins of M 87. 
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They found that the Si-group elements (Si, S, Ar, Ca) and Fe have a stronger peak in the 
central region than O, Ne, and Mg. The SN la enrichment in the centre is characterised by 
a Solar ratio of Si-group elements to Fe. The inferred SN la enrichment in the outer region, 
associated with the ICM of the Virgo cluster, has a lower ratio of Si-gro up elements to Fe 
by 0. 2 dex, which is more cha racteristic for the ICM of other clusters ( iFinoguenov et al.l 
l2000h . lFinoguenov et al.l ( l2002h also find a Ni/Fe ratio of ~1.5 in the central region, while 
they point out that this ratio for other clusters determined by ASCA is ~3. They conclude 
that these abundance patterns confirm the diversity of SN la and to explain the SN la metal 
enrichment in clusters both deflagration and delayed detonation scenarios are required. The 
high SN la yield of Si-group elements in the centre may imply that the Si burning in SN la is 
incomplete. This favours the delayed-detonation models with lower density of deflagration- 
detonation transition, lower C/0 ratio, and lower central ignition density. On the other hand, 
the abundance patterns in the outer region of M 87 are more characteristic of deflagration 
supernova models. 

A bimodal distribution is observed in the magnitude decline of the SN la events. While 
a population of SN la with a slow decline is more common in spiral and irregular galaxies 
with recent star formation, a fain ter and more rapidly decaying populatio n of SN la is more 
common in early-type galaxies jHamuy et al.lll996l : llvanov et al.H 2000Y The SN la with 
incomplete Si burning enriching the central p art of M 87 could be as sociated with optically 
fainter SNe la observed in early- type galaxies. iMannucci et al.l l l2006h found that the present 
data on SN la rates in different populations of parent galaxies indicate that SN la have a 
bimodal delay time distribution. They conclude that about 50 % of SN la explode soon 
after the formation of the progenitor binary system, in a time of the order of 10* years, 
while the delay time distribution of the remaining 50 % of SN la can be described by an 
exponential function with a decay time of ~3 Gyr. The supernovae with short delay times 
might be associated with the supernovae with the slower decline, while the supernovae with 
long delay times might be associated with the rapidly decaying SNe la. If this association is 
true, then the observed Si/Fe of the accumulated SN la ejecta will be a function of the age 
of t he system. 

[Matsushita et al.l ( l2007ah show that SN la ejecta have a higher Fe/Si ratio in regions with 
larger iron-mass-to-light ratio (IMLR). Higher IMLR corresponds to longer accumulation 
time scale. The correlation betwee n the inferr ed Fe/Si rat io of SN la eje cta and the IMLR 
supports the suggestion by Finoguenov et al.l 12002) and iMatsushita et a l. ( 2003) that the 
ave rage Fe/Si ratio produc ed by SN la depends on the age of the system. 

iBohringeretal .1 ( 1200^ compiled the available observational results of the 0/Si/Fe from 
XMM-Newton, Chandra, and ASCA to show that the observed abundance patterns in clus- 
ters are mo re consistent with t he W DD models, which provid e larger Si/Fe ratios than the 
W7 model. iBuote et aT I ( l2003l) and lHumphrev & Buotd ( 12006) also discuss the higher a- 
element enrichment observed in groups and in elliptical galaxies. They conclude that there 
is an increasing evidence that the O is over-predicted by the theoretical SNcc yields, or there 
is a source of ff-element enrichment in addition to that provided by SNcc and SN la. 

Using deep XMM-N ewt on observations of th e clusters 2A 0335-H096 and Sersic 159- 
03. IWemer eFai] (l2006b') and'd e Plaa eTaD (12006^ determined the abundances of 9 elements 
and fitted them with nucleosynthesis models for supernovae. They found that ~30 % of all 
supernovae enriching the ICM were SN la and ~70 % were SNcc- They also found that 
the Ca abundance in these clusters is higher than that expected from the models. Later, 
Ide Pla a et al. ( 2007) used a sample of 22 clusters observed with XMM-Newton to deter- 
mine the abundances of Si, S, Ar, Ca, Fe, and Ni in the ICM within the radius of O.IR^oq 
and they compared the best-fitting abundances with theoretically predicted yields of dif- 
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Fig. 9 Left panel: Comparis on of the [Fe/ Si] ratio of the ICM in the Centauras cluster and in M 87, plotted 
against [Fe/O], from Matsushita et al. 1 2007a). Abundances ar e here on the scale of Feldman 1 1992). The open 
square shows the average value for Galactic metal poor stars iClementini et al. 1999) and the asterisk shows 
the abundance ratios for SNcc calculated by No moto et alj )l997h . The solid line and dashed line show the 
abundance pattern synthesised by SN la with Fe/Si= 1 . 1 (best fit relation for the core of M 87) and Fe/Si=2.6 
(the ratio produced by the W7 model Nomoto et al. 1984), respectively. Right panel: Reconstruction of 
the abundanc e patterns (at the proto-solar Lodders 2003 scale) observed in the sample of 22 clusters of 
Ide Plaa et alj (2007) with theoretical supernova yields. The abundance values indicated by filled circles are 
obtained with XMM-Newton EPIC. The abundances of O and Ne are the average values obtained for Sersic 
159-03 and 2A 0335+096 obtained with R GS. The S N la yields are for a delayed detonation model calibrated 
on the Tycho s upernova remnant by .Badenes et al] (2W6) . The SNcc yield s are for progenitors with Solar 
metallicities bv lNomoto et alj j200d) integrated over the Salpeter IMF. From lde Plaa et al] i2007l) . 



ferent SN la and SNcc m odels (see Fig. [9l). B ecause the SNcc do not have a significant 
impact on the Ar/Ca ratio, Plaa et alj i 2007h used the Ar and Ca abundances in clusters 
as a criterion for determining the quality of SN la models. They found that the Ar and Ca 
abundances in the ICM are inconsistent with currently favored SN la models. However, they 
showed that their best fit Ca/Fe and Ar/Ca ratios are consistent with the empirically modi- 
fied delayed-detonati on SN la model, wh ich is calibrat e d on t he Tycho supernova remnant 
llBadenesetalJl2006i) . Using this model, Ide Plaa et al ] ( l2007l) obtained a number ratio of 
A^SNia/A'sNia+SNcc = 0.44±0.10, which suggcsts that binary systcms in the appropriate mass 
range are very efficient (~ 5 - 16 %) in eventually forming SN la explosions. 

Based on Suzaku observ ations of 2 c lusters (A1060 and AWM 7) and 2 groups of galax- 
ies (HCG 62 and NGC SOT'l. lSato et"aD ( FoOTh ) obtained abundance patterns for O, Mg, Si, 
S, and Fe in the region out to O.Srigo. Contrary to most of the previously mentioned results, 
they found a better fit to the observed abundance patterns using the W7 SN la model rather 
than th e WDD models. Assuming the W7 model for SN la, the yields from Nomoto e t alj 
( l2006h and the Salpeter IMF for SNcc, the ratio of occurrence numbers of SN^; to SN la is 
~3.5. It corresp onds to the nu mber ratio of A'sNia/^SNia+SNcc = 0-22, which is smaller than 
that obtained bv lde Plaa et"^ (12007). The number ratios were determined only considering 
the metals in the ICM and metals locked up in the stars were not considered. 

6.2 Summary of the efforts of using ICM abundances to constrain supernova models 

The abundance patterns obtained during many deep observations of clus ters of galaxies 
with XMM-Newton clearly favor the delayed-detonation SN la models (see lBohringer et aU 
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l2005h and the observed abundance ratios of Ar/Ca are relatively well reproduced using 
empirically modified delayed-detonation models calibrated on the Tycho supernova rem- 
narit (de Plaa et al. 2007). However, several observations still favor the W7 SN la model 
dSato et aI.II2007b ) and in most of the cases the abundanc es are not determ i ned \y ell enough 
to really discriminate between the diff^erent models. As Finog uenov et al.l J2002h suggest It 
might even be possible that both the deflagration and delayed-detonation types of SN la play 
a role in the chemical enrichment of clusters of galaxies. 

XMM-Newton RGS observations of the hot gas in clusters of galaxies and in elliptical 
galaxies reveal Ne/0 ratios which are higher than Solar jXu et al .l2002lJPeterson et al.l2003l : 
IWemer et al.l2006al lbh. SN la produce very little O and Ne, and their ratio can be used to put 
constraints on the progenitors of SNcc- From the SNcc models discussed in Sect.[2]only the 
models with pre-enriched progenitors predict higher than Solar Ne/0 ratio. 

While it is well known that elements from O up to the Fe-group are primarily produced 
in supernovae, the main sources of C and N are still being debated (see Sect. Deep 
observations of nearby X-ray bright elliptical galaxies and bright clusters with the RGS on 
XMM-Newton allow to measure the C and N abundances in the hot gas. By combining the 
RGS data obtained in two deep XMM-Newton observations of M 87. IWemer et al.l l l2006al) 
determined relatively accurate C and N abundance values. The small 0/Fe ratio and large 
C/Fe and N/Fe ratios found in M 87 suggest that the main sources of C and N are not the 
massive stars that also produce large quantities of O, but the low- and intermediate-mass 
asymptotic giant branch stars. 



7 Abundances as function of cluster mass and redshift 

lAmaud et all ( Il992h studied the correlations between the gas mass, Fe mass, and optical 
luminosity in clusters of galaxies. They found that the gas mass in the ICM is highly corre- 
lated and the Fe mass is directly proportional to the inegrated optical luminosity of elliptical 
and lenticular galaxies. They found no correlation with the integrated optical luminosity of 
spiral galaxies. The ratio of the gas mass in the ICM to the stellar mass increases with the 
clu ster richness. 

lRenzinil ( ll997h showed that the IMLR is very similar in all clusters of galaxies with ICM 
temperature above ~1 keV. However, below ~1 keV the IMLR seems to drop by almost 3 
orders of magnitude. This drop at low temperatures is due to the combination of lower gas 
mass and lower measured Fe abundance in cool systems. However, it was shown that the low 
Fe abundances measured in many groups were artifacts of wrong temperature modelling (see 
late r). 

iRenzinil l ll997l , l2004h showed that clusters hotter than ~2.5 keV have similar abun- 
dances, with very small dispersion. Below 2.5 keV, the best fit metallicities show a large 
range of values. The previously reported very low abundance values for low temperature 
groups and elliptical galaxies turn out to be much higher when two- or multi-temperature 
fitting is applied (see Sec t. 14.11 and iBuote & FabiarJI 19981 : lBuotell2000bl : | Buote et al. I I2OO3I : 
iHumphrev & Buote" 2006V However, the outer regions of galaxy groups still show very low 
metallicities. An XMM-Newton observation of the group NGC 5044 indicates an iron abun- 
dance in the hot gas of Z pe ~ 0.15 Solar in the region between 0.2 and 0.4 of the virial 
radius jBuote etaUl 20041) . This result shows that the total iron mass to optical light ratio in 
NGC 5044 is about 3 times lower than that observed in rich clusters. 

We compile d the abundance valu es measured within 0.2/?5oo for the sample of 22 clus- 
ters analysed by Ide Plaa et al.l | |2007|) . In the top panel of Fig. [TO] we show the best fit Fe 
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Fig. 10 Top panel: Best fit Fe abundance as a f unctio n of ICM temperature determined witliin 0.2i?5oo for a 
sample of clusters investigated by de Plaa et"ar| j2007|). We compare the observed Fe abundances with values 
from simulations i Borgani et al. 2008 - Chapter 18, this volume) for a Salpeter IMF, extracted within the 
same physical radius. Bottom panel: Abund ance ratios of Si , S, Ca, and Ni with respect to Fe as a function of 
the cluster temperature, compiled from ide Plaa et alj j2007ft . The abundance ratios are consistent with being 
constant as a function of cluster mass. Cooling core clusters are shown as red circles and non-cooling core 
clusters as blue tr iangles. Abundances in both panels are plotted with respect to the proto-solar abundance 
values of lLoddersI r2003i) . 
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Fig. 11 Left panel: Mean Fe abundance from combined fits within six redshift bins (red circles) and the 
weighted average of single source observations in the same bins (black square) from Balestra et al. 1 2 0071) . 
The shaded areas show the rms dispersion. The dashed line shows the best power-law fit over the redshift 
bins. Right panel: Mean ICM metal abundances in eight redshift bins within R^qq with and without the 
ce ntral O . I5R50 excluded (from lMaughan et alj200l) . The solid line shows the sup ernova emichment model 
of lEttoril <2005h . Both panels use the abundance scale of lAnders & Grevessd jl989l) . 



abundance as a function of ICM temperature. Our data confirm earlier results showing that 
while in cooler clusters, in the temperature range of 2-A keV, Zpe shows a range of values 
between 0.2-0.9 Solar, in hot massive clusters (kT > 5 keV) the Fe abundance is on average 
lower and equal to Zpe ~ 0.3 Solar. We compare the observe d trend with values from sim- 
ulations extracted within O.lRsoo as the XMM-Newton data (Borg ani et al ] l2008l - Chapter 
18, this volume). The data are simulated for the Salpeter IMF. The data are more consistent 
with the simulated values at the lower temperatures. At higher temperatures the simulated 
abundances are higher than the observed values. The observed trend might be linked to a 
decrease of the star formation efficiency with increasing cluster mass (Lin et al. 2003). 

In the bottom panel of Fig.[TO]we show the abundance ratios of Si, S, Ca, and Ni wit h 
respect to Fe as a function of the cluster temperature, compiled from lde Plaa et al. 112003). 
The abundance ratios are consistent with being constant as a function of temperature, which 
is related to the cluster mass. The intrinsic spread in abundance ratios between clusters i s 
smaller than 30 %. This finding contradicts the ASCA results ( iBaumgartner et al" 
which show an increase of the Si/Fe ratio from 0.7 to 3 Solar in the 2-8 keV temperature 
range. This result indicates that the ratios of SN la and SNcc contributing to the enrichment 
of the ICM are similar for clusters of all masses. 

Investigating a sample of 56 clusters at z > 0. 3. iBalestra et"aD ( l2007h confirmed the 
trend of the Fe abundance with cluster temperature also in the higher redshift clusters. They 
show that the abundance measured within (0.15-0.3)/?vir in clusters below 5 keV is, on 
average, a factor of ~2 larger than in hot clusters, following a relation Z{T) - 0.887'"° '*^. 

iBalestra et all ilOOf) also found an evolution of the Fe content in the ICM with redshift 
(see the left panel of Fig.lllb. They conclude that the average Fe content of the ICM at the 
present epoch is a factor of ~2 larger than at z - 1.2. While at the redshift z > 0.5 they 
observe an average ICM abundance of Zpe - 0.25 Zq, in the redshift range of z - 0.3-0.5 
they measure Zpe - 0.4 Zg. They parametrise the decrease of metallicity with redshift with a 
power-law ~ (1 + z)~' '^^. The evolution in the Fe abundance with redshift was confirmed by 
iMaughan et al]( l2007h , who investigated a sample of 116 clusters of galaxies at 0.1 <z< 1.3 
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in the Chandra archive (see the right panel of Fig. II 111. They found a significant evolution 
with the abundance dropping by 50 % between z ~ 0.1 and z ~ 1- They verified that the 
evolution is still present, but less significant, when the cluster cores are excluded from the 
abundance measurement (they use apertures of 0-1 R500 and 0.15-1 R500), indicating that 
the evolution is not solely due to d i sappea rance of relaxed, cool core clusters at z > 0.5. 

Simulations bv lKapferer et alj ( l2007bh show that between redshifts of z = 1 and z = a 
galaxy cluster accretes fresh gas, with the gas mass increasing during this time period by a 
factor of ^3. Which means that in order to get an increase of the metallicity during this time 
period the freshly accreted gas must be pre-enriched or the input of metals into the ICM in 
this time interval must be larger than expected. 



8 Future of ICM abundance studies 

A major improvement in X-ray spectroscopy will be the employment of micro-calorimeters, 
and especially of the so-called Transition-E dge Sensors (TES) on future satellites (for more 
information on future instrumentation see IPaerels et al.l [20081 - Chaper 19, this volume). 
Compared to the current state-of-the-art CCD technology, this new type of detectors promises 
a factor of > 20 improvement in spectral resolution (< 2 - 5 eV, depending on the optimal- 
isation of the detector). To reach this high spectral resolution, the TES detectors use the 
transition edge between normal conductivity and superconductivity that is present in a num- 
ber of materials. At temperatures of ~0.1 K, the conductivity of the absorbing material is 
very sensitive to the temperature. When an X-ray photon hits the detector, its temperature 
rises and the resistance of the material increases. By measuring the increase of the resis- 
tance, one can calculate with high accuracy the energy of the X-ray photon. The main elTort 
today is to build an array of TES detectors, with imaging capabilities. Next major proposed 
missions with TES arrays on board are, amongst others, the X-ray Evolving Universe Ex- 
plorer (XEUS), Explorer of Diffuse emission and Gamma-ray burst Explosions (EDGE), 
and Constellation-X. 

The combination of the large eff^ective area (~5 m^) and high spectral resolution (< 
2-5 eV) of XEUS will open up new possibilities in cluster abundance studies, with enor- 
mous improvement in the abundance determinations. In the left panel of Fig. [12] we show a 
simulated spectrum of the cluster Sersic 159-03 (from de Plaa 2007), for the planned TES 
detector that will be part of the Narrow-Field Imager (NFI) on XEUS. For comparison also 
the model at the resolution of XMM-Newton EPIC is shown. With XEUS the abundances 
will be measured with an accuracy of ~10"-' Solar, which is more than an order of magnitude 
improvement compared to EPIC. This accuracy might be sufficient to detect the contribu- 
tion of the Population-Ill stars to the enrichment of the ICM. Also spectral lines from more 
elements will be resolved, which will enable better tests for supernova models. 

The large collecting area of XEUS in combination with the CCD-type detectors of the 
Wide-Field Imager will be excellent to observe the spatial distribution of metals, which 
will help us to constr ain the enrichment history of clusters. Simulated metallicity maps by 
iKapferer et all ( l2006h for Chandra, XMM-Newton, and XEUS are shown in the right panel 

of Fig. [m 

The proposed EDGE satellite would also open up new possibilities in the chemical abun- 
dance studies. The combination of the planned long exposure times (of the order of few Ms), 
large field of view (1.4 degrees) and good angular resolution (Half Power Diameter of 15") 
of its CCD type Wide-Field Imager, and the high spectral resolution (3 eV at 0.6 keV) of the 
Wide-field Spectrometer promises new possibilities in the studies of metal abundances in 
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Fig. 12 Top panel: Simulated line spectrum of the inner 30" of the cluster of galaxies Sersic 159-03 for a 
100 ks observation with the TES array planned for XEUS. The continuum model is subtracted. F or com- 
parison the EPIC spectrum of the same source is overplotted as the thick line. From[de Plaa' f2007|). Bottom 
panel: Simulated X-ray metallicity maps as seen with Chandra, XMM-Ne wton, and XE US, respectively. The 
plot in the bottom right comer shows the original simulation results. From lKapferer et al., ( .2006) . 
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clusters. EDGE would not only increase the number of detected spectral lines and improve 
the accuracy of the abundance measurements, but it would also enable us to accurately mea- 
sure the chemical abundances and the spatial distribution of metals in the faint outskirts of 
clusters out to their virial radius, which is beyond the possibilities of current instruments. 
This would provide us with insights of crucial importance for the understanding of the chem- 
ical enrichment processes, nucleosynthesis, and of the metal budget of the Universe. 
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